Search for: All records

Abstract Silicate vapors play a key role in planetary evolution, especially dominating early stages of rocky planet formation through outgassed magma ocean atmospheres. Our open-source thermodynamic modeling software “VapoRock” combines the MELTS liquid model with gas-species properties from multiple thermochemistry tables. VapoRock calculates the partial pressures of 34 gaseous species in equilibrium with magmatic liquid in the system Si–Mg–Fe–Al–Ca–Na–K–Ti–Cr–O at desired temperatures and oxygen fugacities (fO₂, or partial pressure of O₂). Comparison with experiments shows that pressures and melt-oxide activities (which vary over many orders of magnitude) are reproduced to within a factor of ∼3, consistent with measurement uncertainties. We also benchmark the model against a wide selection of igneous rock compositions including bulk silicate Earth, predicting elemental vapor abundances that are comparable to (Na, Ca, and Al) or more realistic than (K, Si, Mg, Fe, and Ti) those of the closed-source MAGMA code (with maximum deviations by factors of 10–300 for K and Si). Vapor abundances depend critically on the activities of liquid components. The MELTS model underpinning VapoRock was calibrated and extensively tested on natural igneous liquids. In contrast, MAGMA’s liquid model assumes ideal mixtures of a limited set of chemically simplified pseudospecies, which only roughly approximates the nonideal compositional interactions typical of many-component natural silicate melts. Finally, we explore how relative abundances of SiO and SiO₂provide a spectroscopically measurable proxy for oxygen fugacity in devolatilized exoplanetary atmospheres, potentially constrainingfO₂in outgassed exoplanetary mantles. 

Abstract MESSENGER observations suggest a magma ocean formed on proto-Mercury, during which evaporation of metals and outgassing of C- and H-bearing volatiles produced an early atmosphere. Atmospheric escape subsequently occurred by plasma heating, photoevaporation, Jeans escape, and photoionization. To quantify atmospheric loss, we combine constraints on the lifetime of surficial melt, melt composition, and atmospheric composition. Consideration of two initial Mercury sizes and four magma ocean compositions determines the atmospheric speciation at a given surface temperature. A coupled interior–atmosphere model determines the cooling rate and therefore the lifetime of surficial melt. Combining the melt lifetime and escape flux calculations provides estimates for the total mass loss from early Mercury. Loss rates by Jeans escape are negligible. Plasma heating and photoionization are limited by homopause diffusion rates of ∼10⁶kg s⁻¹. Loss by photoevaporation depends on the timing of Mercury formation and assumed heating efficiency and ranges from ∼10^6.6to ∼10^9.6kg s⁻¹. The material for photoevaporation is sourced from below the homopause and is therefore energy limited rather than diffusion limited. The timescale for efficient interior–atmosphere chemical exchange is less than 10,000 yr. Therefore, escape processes only account for an equivalent loss of less than 2.3 km of crust (0.3% of Mercury’s mass). Accordingly, ≤0.02% of the total mass of H₂O and Na is lost. Therefore, cumulative loss cannot significantly modify Mercury’s bulk mantle composition during the magma ocean stage. Mercury’s high core:mantle ratio and volatile-rich surface may instead reflect chemical variations in its building blocks resulting from its solar-proximal accretion environment.